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NOMENCLATURE 


A  Helmholtz  free  energy 
C  number  of  components 

F  Gibbs  free  energy 

g  acceleration  due  to  gravity 

h  height  above  a  free  surface 

H  enthalpy 

L  latent  heat 

Lt  latent  heat  of  evaporation 

m  mass 

Nv  Nb  mole  fractions 
p  pressure 

p\,  Pj  pressure  of  phases  i  and  j 

pQ  reference  pressure 

P  number  of  phases 

Q  heat 

r1 ,  r2  radii  of  curvature  of  two  mutually  perpen¬ 
dicular  sections 

rp  radius  of  particle 

rv  radius  of  a  liquid  vein 

rfHm  mean  radius  of  curvature  of  film 

Tjj  mean  radius  of  curvature  of  phase  boundary 

/?1  ,R2  principal  radii  of  curvature 

5  entropy 

5j  entropy  of  /th  phase 


T  temperature  (K) 

/j  temperature  of  /th  phase  (K) 

Tm  meeting  temperature  (°C) 

f0  reference  temperature  (K) 

U  internal  energy 

y  specific  volume 

specific  volume  of  /th  phase 
V  volume:  also  molar  volume 

W  work 

xa  mole  fraction  of  solvent 

xb  mole  fraction  of  solute 

6  contact  angle 

p  chemical  potential 

Pj  chemical  potential  of /th  phase 

it  osmotic  pressure 

p  density 

pj  density  of /th  phase 

Ojj  interfacial  energy  between  /th  and/th  phases 
<f>  solid  angle  of  portion  of  surface 

Subscripts 

Q,bfiJ  indices  for  phases  or  components 

g,2,s  gas,  liquid,  solid,  respectively 
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INTRODUCTION  TO  THE  BASIC 
THERMODYNAMICS  OF  COLD 
CAPILLARY  SYSTEMS 

Samuel  C.  Colbeck 


INTRODUCTION 

The  concepts  of  thermodynamics  apply  to  a  wide 
variety  of  cold  regions  materials  in  which  radius  of 
curvature  effects  are  important  (for  example,  frozen 
soils,  snow,  and  ice  fog).  The  basic  principles  applic¬ 
able  to  these  and  many  other  important  phenomena 
are  reviewed  here,  and  some  example  applications  are 
given. 

This  report  is  intended  to  serve  as  a  brief  introduc¬ 
tion  or  review  of  the  elements  of  thermodynamics 
which  are  applicable  to  cold  regions  materials.  Refer¬ 
ences  to  additional  sources  of  information  are  pro¬ 
vided,  since  only  an  outline  of  thermodynamics  is 
given  here. 

THERMODYNAMIC  PRINCIPLES 

J.  Willard  Gibbs  (» 839-1 903)  greatly  increased  our 
understanding  of  thermodynamics,  especially  the 
thermodynamics  of  reversible  processes  (Gibbs  1928). 
Some  of  his  most  important  contributions  are  re¬ 
viewed  here. 

The  first  law 

The  first  law  of  thermodynamics  expresses  energy 
balance  as 

dU=dQ+dW  (1) 

where  dU  is  the  change  in  internal  energy  of  the  sys¬ 
tem,  dQ  is  a  small  amount  of  heat  energy  entering  the 
system,  and  dW  (usually  taken  as  -pdV)  is  the  work 
done  on  the  system  by  a  change  of  volume  dV against 

the  ores* urep,  This  equation  is  inexact  since  f*pdV 


depends  on  the  path  followed  in  going  from  A  to  B 
(except  in  reversible  processes,  where  the  integral  is 
independent  of  path).  We  define  enthalpy  H  as 

H  =  U+pV.  (2) 

Hence,  another  form  of  the  first  law  is 

bQ^dH-Vdp.  (3) 

The  second  law 

From  the  second  law  of  thermodynamics 

f  42  <  0  (4) 

7  T 

where  the  integration  is  over  a  closed  path  (that  is, 
the  integration  starts  and  stops  at  the  same  point). 

The  inequality  applies  if  the  process  is  irreversible  (for 
example,  losses  due  to  friction  or  strain  energy)  and 
the  equality  applies  if  the  process  is  completely  revers¬ 
ible  (for  example,  the  adiabatic  expansion  and  contrac¬ 
tion  of  a  perfect  balloon).  When  a  reversible  process 
occurs,  the  integral  must  be  independent  of  the  path 
followed  from  A  to  B.  We  define  entropy  5  as 

VSA  >  fA*  (S) 

where  A  and  B  are  any  two  positions.  Again,  the  in¬ 
equality  applies  to  an  irreversible  process.  Therefore, 
entropy  must  always  increase  in  any  irreversible  process 
in  an  isolated  system.  (In  other  words,  the  world  is 
filling  up  with  entropy!) 


dU+pdV-TdS  =  0 


(15) 


Entropy  maximum 

Unless  the  system  is  in  equilibrium,  entropy  will 
be  increasing  according  to  eq  5.  Therefore,  entropy 
is  at  its  maximum  at  equilibrium.  From  eq  1  and  5, 

55  >  b~  P-h  y  (6) 

T 


or 


t>U+phV-TbS  <  0. 


(7) 


we  find  a  very  general  relation 

dF  -  Vdp-SdT.  (16) 

In  a  single-component  system  of  two  or  three 
phases  (for  example,  vapor,  liquid  and  solid  H20), 
eq  16  can  be  applied  to  each  phase  or  to  the  whole 
system  at  once.  When  the  phases  are  in  equilibrium, 
the  free  energies  of  the  phases  are  equal. 


At  constant  volume  (6  V  ~  0), 

8U-TbS<  0  (8) 

or  at  constant  volume  and  temperature, 


Chemical  potential 

If  more  than  one  component  exists  (for  example, 
H20  and  Si02),  eq  16  can  be  generalized.  Gibbs  de¬ 
fined  the  chemical  potential  (or  partial  molal  free 
energy)  p  as 


6(U-TS)<0.  (9) 

The  expression  U-TS  is  the  definition  of  the  Helm¬ 

holtz  free  energy  A ,  so  8A  <  0.  Accordingly,  a  sys¬ 
tem  will  move  toward  lower  values  of  A,  and  at  equi¬ 
librium  A  will  assume  a  minimum. 

If  the  system  is  at  constant  temperature  and  pres¬ 
sure  but  variable  volume,  then  eq  7  is  rewritten  as 

b(U+pV-TS)  <  0.  (10) 

The  Gibbs  free  energy  F  is  defined  as 

F=U+pV-TS  (11) 


or 


F=H-TS  (12) 

or 

F  =  A+pl/.  (13) 

In  this  system  (6p  =  5T=  0),  the  Gibbs  free  energy  is 
at  its  minimum  at  equilibrium.  It  follows  that  when 
liquid,  gaseous,  and/or  solid  phases  are  in  equilibrium, 
their  free  energies  are  equal.  If  they  are  not  equal, 
the  phase  with  the  largest  value  of  F  will  give  up  mass 
to  the  other  phases  until  equilibrium  is  established. 

Maxwell's  relations 

We  can  derive  a  very  useful  form  of  the  Gibbs  equa¬ 
tion.  Differentiating  eq  1 1 


Mab  = 


9wab 


(17) 


where  m  is  the  number  of  moles  per  unit  volume,  a  is 
the  component,  b  is  the  phase.  Phase  changes  occur 
spontaneously  such  that  material  is  transferred  from 
a  phase  at  a  high  potential  to  a  phase  at  a  low  poten¬ 
tial.  Equilibrium  occurs  when  the  potential  of  each 
component  is  the  same  in  each  phase.  In  general, 


dF  =  ~SdT+  Vdp+Xp^dm^  (18) 

Gibbs-Duhem  equation 

In  another  form,  called  the  Gibbs-Duhem  equation 
(i dF =  0  at  equilibrium), 


SdT-Vdp+  =  0  .  (19) 

In  a  system  of  ice,  liquid  water  and  moist  air,  there  are 
two  components  (air  and  water),  but  the  air  component 
can  be  ignored.  (If  it  is  included,  only  the  partial  pres¬ 
sure  of  the  dry  air  need  be  specified.  Air  dissolved  in 
the  liquid  phase  decreases  the  melting  temperature  by 
0.0024°C  when  the  liquid  is  saturated.)  Then  the 
Gibbs-Duhem  equation  reduces  to  a  simple  relation 
which  is  commonly  used, 

dp  ~  vdp-SdT  (20) 


where  pt  v  and  5  are  expressed  per  unit  mass.  Equa¬ 
tion  20  can  be  applied  to  each  phase  separately.  Then 

dp ,  =  v{dprS^dT  (21) 


dF  =  dU+pdV+  Vdp-TdS-S  dT  (14) 

and  rewriting  eq  7  for  a  reversible  process 


where  /  is  the  solid,  liquid  or  vapor  phase.  When  only 
one  component  exists, 
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(22) 


(24) 


dml 


At  equilibrium,  F  is  at  its  minimum,  so  mass  transfer 
will  occur  between  the  nhases  (mass  will  transfer 
away  from  phases  of  high  potential)  until  the  poten¬ 
tial  is  the  same  for  all  phases.  Thus,  at  equilibrium 
between  phases  /  and  /, 

(23) 

Phase  rule 

Gibbs’  phase  rule  tells  us  how  many  variables  must 
be  specified  at  equilibrium.  When  there  are  no  chemi¬ 
cal  reactions,  and  C  components  (for  example  water, 
air,  Si02,  salts)  and  P  phases  (for  example,  solid, 
liquid,  vapor)  are  considered  independently,  there  are 
/’(C+l)  variables  which  must  be  specified  at  equilibri¬ 
um.  Since  the  phases  are  in  equilibrium,  there  are 
C(P- 1 )  restrictions  because  the  chemical  potentials 
must  be  equal.  Also,  since  the  temperatures  of  each 
phase  must  be  equal,  there  are  P- 1  additional  restric¬ 
tions.  Thus,  for  a  system  with  curved  interfaces, 
there  are  P(C+ 1  )-C(/*- 1 )  -  (P- 1 )  or  C+ 1  variables 
which  must  be  specified  at  equilibrium.  This  is  the 
number  of  “degrees  of  freedom”  or  the  “variance”  of 
the  system.  Note  that  for  flat  interfaces  there  are 
P-1  additional  restrictions,  since  the  pressures  must 
also  be  equal.  Thus,  for  non-capillary  systems,  there 
are  C+ 2 -^degrees  of  freedom. 


where  R}  and  R2  are  the  principal  radii  of  curvature  of 
the  surface  separating  these  two  phases.  The  mean 
radius  of  curvature  (rjj)  is  defined  as 


(25) 


Consider  an  element  of  length  of  the  perimeter 
at  A .  The  component  along  PN  of  the  force  on  the  ele¬ 
ment  at/4  is  ajj5C  sin 0  where  a}j  is  the  surface  energy 
between  phases  /  and  j.  Since  0  is  small, 


sin 0  ~  0 


(26) 


Thus,  the  force  along  PN  due  to  this  element  is 
Ojj  (p/r1  )5B  and  is  positive  in  the  direction  of  phase  /. 
At  points  C  and  D  there  are  similar  forces  equal  to 
°\]  (p/^)^-  Integrating  around  the  perimeter  of 
this  section,  we  find  that  the  total  force  directed  to¬ 
wards  phase  /  is  jrp2  [(1 /rj)+ (l/r2)]  •  This  must 
be  balanced  by  a  pressure  difference  between  the 
phases,  or 


SURFACE  FORCES 

Since  capillary  systems  have  curved  interfaces  be¬ 
tween  the  solid,  vapor  and/or  liquid  phases,  pressure 
differences  can  occur  among  the  phases.  These  pres¬ 
sure  differences  must  be  included  in  the  derivation  of 
the  governing  equations.  Laplace’s  equation  describes 
the  pressure  difference  between  two  adjoining  phases; 
it  is  derived  here  because  of  its  basic  importance  in 
cold  capillary  systems.  The  pressure  difference  arises 
from  the  tension  in  the  interfacial  layer  (surface  free 
energy  or  surface  tension)  and  from  the  curvature  of 
the  interface  in  much  the  same  way  that  a  stretched 
balloon  supports  a  pressure  difference  between  the 
inside  air  and  the  outside  air. 

Figure  1  depicts  a  curved  surface  between  two 
phases  (/  and  /).  Let  the  distance  along  the  surface 
from  the  point  P  to  the  edge  of  the  surface  be  every¬ 
where  equal  to  p.  Let  the  two  sections  A  BN  and  CDN 
be  mutually  perpendicular.  Then  the  radii  of  curva¬ 
ture  in  these  two  sections  (^  and  r2)  completely  de¬ 
scribe  the  surface  separating  the  phases  /  and  /.  Ac¬ 
cording  to  Euler’s  theorem, 


Figure  h  Portion  of  the  surface  dividing  phase  i  (be¬ 
low)  from  phase  j  (above). 
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Vapor 


Figure  2.  Force  equilibrium  at  a  point  of  liquid-vapor  contact 
on  a  solid  surface. 


or 


and 


p~P) 


or  finally,  Laplace’s  equation, 


Pi-Pi 


(28) 


(29) 


where  is  the  an  radius  of  curvature  between 
phases  /  and  j.  By  convention,  we  call  positive 
when  phase  /  is  on  the  concave  side  of  the  interface. 

Now  consider  the  line  of  contact  of  three  adjacent 
phases.  Summing  forces  along  the  tangent  to  the  sur¬ 
face  at  the  point  of  contact  P  in  Figure  2, 


dpq  =  vv  dp^-S^dT^. 

The  dp* s  vanish  at  equilibrium.  For  a  flat  surface, 
pc  =ps,  and  from  eq  5,  the  difference  in  entropy  at 
constant  temperature  is  related  to  the  latent  heat  of 
fusion  L  (usually  assumed  to  be  constant)  by 


5C-Ss  =  LIT. 


03) 


Thus, 


(vs-vt)dp+{Si-S%)dT=0 


or 


asg  =  ase+agecos(?-  <3°) 

This  is  Young’s  equation.  Since  for  most  materials 
the  contact  angle  (0)  is  small  (that  is,  cos0  -  1 ), 


(31) 


or  £sg  is  generally  much  greater  than  either  osC  or  agC. 
Surface  energies  osg  and  osC  cannot  be  measured  di¬ 
rectly  but  must  be  deduced  by  methods  such  as  meas¬ 
urement  of  6  and  ogC  or  grain  boundary  observation. 

There  is  some  controversy  over  the  correct  inter¬ 
pretation  of  osg  and  asC.  In  the  derivation  of  Laplace’s 
equation  (eq  29),  we  assumed  that  only  isotropic 
stresses  (or  pressures)  occur  in  the  phases.  Therefore, 
the  equation  may  only  be  applied  in  the  absence  of 
deviatoric  stress;  application  to  solids  must  be  done 
carefully. 


GENERAL  APPLICATIONS 

Clausius-Clapeyron  equation 

For  a  liquid  and  solid  with  a  flat  surface  in  com¬ 
mon,  eq  21  shows  that 


(v^-v^dp-  L(dT/T). 

Integrating  from  p  =  0  (gauge)  and  T  -  T0  (273.1 7K) 
and  assuming  that  L  is  not  dependent  on  pressure  or 
temperature, 


where  T  and  T0  are  in  kelvjns.  We  define  the  melting 
temperature  Tm  in  °C  as 

7m(°C)  =  r(K)-70(K).  (35) 


Now 


,n  —  =  In  hslR 
To  To 

*Mi  +  rm/T0) 

or  since  7*m  «T0 , 

lnZ  =  Z"L. 

To  T0 


(36) 


dPi  ~  vidpi-SidTi 


(32) 


Combining  eq  34  and  36  to  get  the  Clausius-Clapeyron 
equation, 
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■iipayw 


5>  (i  -  ±), 

L  Ve  Pj 


(37) 


ln  Pg  =  .  v±  2qg« 
P0  *7"  rge 


(43) 


where  the  density  p  is  the  inverse  of  the  specific  vol¬ 
ume  i/.  For  ice  in  contact  with  liquid  water, 

rm  -  -0.00736p  (°C/bar). 

Since  pc  >  ps,  the  melting  temperature  of  water  de¬ 
creases  with  increasing  pressure.  This  is  probably  the 
best  known  result  of  phase  equilibrium  thermody¬ 
namics. 


In  a  capillary  system,  smaller  pores  have  smaller  (but 
positive)  values  of  rgC  and  lower  vapor  pressures.  Thus 
the  larger  pores  are  preferential  sites  for  evaporation. 

Capillary  rise 

In  a  capillary  system,  gravity  must  be  considered, 
because  the  system  height  is  large  compared  to  the 
size  of  a  meniscus.  When  the  liquid  is  connected 
throughout  the  pore  space, 


Kelvin’s  equation 

For  liquid  and  vapor  in  contact  at  a  constant  tem¬ 
perature  at  equilibrium,  eq  21  gives 

Vgdp^^dpz 

v%dpg  =  v^dp^-v^dp^v^dp  g  (38) 

(i/g  -vi)dpg  =  vz(dpQ-dpg). 

Using  vg  »  v c  and  Laplace’s  equation  (eq  29), 


where  r%g  is  positive  when  measured  from  the  liquid 
side  of  the  interface. 

For  one  mole  of  a  perfect  gas, 

Pg  =  RT  (40) 

so 


Pa~Po-Ps9h  f44) 

where  h  is  the  height  above  a  free  surface  (p0  =  pc  =pg 
at  the  free  surface)  and  g  is  the  acceleration  due  to 
gravity.  From  eq  21  for  a  capillary  system  at  a  con¬ 
stant  temperature  at  equilibrium, 

vgdpg  =  vQdpQ 

so 

vgdpg  =  v9(-psgdh). 

Using  =  (1/p8)  and  eq  40, 

RT^S.  =  - gdh . 
p% 

Integrating  from  0  to  h  and  p0  to pg, 

In^i  =--£-/»  (45) 

Po  PT 


For  a  flat  interface,  /-fig  is  infinite  and  the  pressure  is 
p0.  Thus  eq  41  can  be  integrated  to  give  Kelvin’s 
equation 

in  P-±  =  H  .  (42) 

Po  *T  r9 g 

This  shows  that  the  vapor  pressure  increases  asrCg  de¬ 
creases.  Thus,  large  droplets  of  water  in  the  atmo¬ 
sphere  have  lower  vapor  pressures  than  small  droplets, 
and  the  larger  droplets  tend  to  grow  at  the  expense  of 
the  smaller  droplets  (that  is,  the  small  droplets  evapo¬ 
rate  because  of  their  large  vapor  pressure,  while  the 
vapor  condenses  on  the  large  droplets). 

If  the  radius  of  curvature  r  is  measured  from  the 
air  side  of  the  interface, 


or 

pg^  p0exp(-gh/RT).  (46) 

Thus  the  vapor  pressure  decreases  with  height  above  a 
free  surface.  This  is  an  alternate  form  of  Kelvin's 
equation. 

Consider  the  vapor  pressure  pg  as  constant  in  a 
capillary  system  and  letpg  =pB  =p0  at  a  free  surface. 
Above  the  free  surface,  pc  decreases  according  to  eq 
44,  while pg  is  nearly  constant.  Then 

Pg-P«  =Pg-(P0-Pfi^) 

or 

p*-p« =Pt^- 

Since  p9gh  >  0,pg-ps  >  0.  By  Laplace’s  equation, 
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2 agc /rgc  '  hence,  rgC  >  0.  In  other  words,  the 

vapor/liquid  interface  must  have  a  positive  radius  of 


Menisci  and  grains 

In  a  meniscus  between  two  rounded  particles  (Fig. 
3)>  rge  >  °f  but  a  liquid  film  cannot  cover  the  particle 
without  violating  this  condition.  The  mean  curvature 
is  given  by 

2-11  mo\ 


where  and  r2  are  the  principal  radii.  Both  are  posi¬ 
tive  values  as  long  as  1  /r  2  is  subtracted  from  1  /r1 . 
Remember,  /gC  is  measured  from  the  vapor  side. 
However,  if  a  film  of  water  surrounds  the  grain  of  ra¬ 
dius  rp,  the  mean  radius  of  the  film  is  given  approxi¬ 
mately  by 


rP  rP 


'film  =  -  V  <5°) 

Thusrfj|m  <  0,  which  violates  the  condition  >  0. 
Therefore,  if  a  film  exists,  it  must  not  have  the  proper¬ 
ties  of  a  liquid;  that  is,  it  can  support  elastic  stresses 

Dissolved  impurities 

When  another  component  exists  in  the  form  of  a 
dissolved  solid,  eq  19  must  include  its  mass  and  chemi¬ 
cal  potential.  At  equilibrium  the  chemical  potentials 
of  the  two  components  are  equal,  but  the  mass  of  the 
dissolved  solid  must  still  be  considered.  In  effect, 
another  term  is  added  to  equations  like  eq  37  and  42. 
This  term  reduces  the  equilibrium  temperature  by  an 
amount  which  depends  on  the  quantity  of  the  impur¬ 
ity  and  the  particular  geometry  of  the  system. 

In  soil  science  it  is  popular  to  describe  the  effect  of 
dissolved  salts  using  the  term  “osmotic  pressure/’  The 
effect  of  salts  (at  constant  temperature  and  pressure) 
is  to  decrease  the  free  energy.  This  decrease  can  be 
restored  at  a  constant  temperature  by  raising  the  pres¬ 
sure  an  amount  tt;  that  is,  a  pure  liquid  will  be  in  equi¬ 
librium  with  a  solution  if  the  solution  is  at  an  osmotic 
pressure  n  higher  than  that  of  the  pure  liquid. 

For  a  pure  liquid,  the  free  energy  change  with  pres- 


AF  =  RT  In pc/p0.  (51) 

In  general,  the  dependence  of  free  energy  at  constant 


Figure  3.  Liquid  meniscus  between 
two  rounded  particles. 


temperature  is  expressed  by 


where  V  is  the  volume  per  mole.  Thus  for  a  pressure 
increase  tt, 

AF  =  Vn.  (53) 

Adding  eq  51  and  eq  53  so  that  the  decrease  (due  to 
salinity)  and  the  increase  (due  to  pressure)  balance, 

”  +  ^Hp/po)  =  0-  (54) 


From  Raou/t’s  observations, 
P/Po  = 


where  xa  is  the  mole  fraction  of  solvent.  The  total 
fraction  is 


1  =*a  +  *b 


where  xb  is  the  mole  fraction  of  solute.  For  dilute 
solutions  it  immediately  follows  that  the  ratio  of  ihe 
mole  fractions  N  is 

lnxaft  -Nb/Na  (57) 

so  7r  =  RTN^KN^  V)  or  the  osmotic  pressure  is  propor¬ 
tional  to  the  number  of  moles  of  solute  per  unit  vol¬ 
ume  of  dilute  solution  at  a  constant  temperature. 


COLD  REGIONS  APPLICATIONS 
Frozen  soils 

A  number  of  important  topics  have  already  been 
discussed,  including  capillary  condensation,  capillary 
rise,  meniscus  shape  and  ionic  impurities.  In  these 
cases  we  assumed  constant  temperature  and  examined 
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Figure  4.  Coexisting  ice -  and  liquid-filled  Figure  5.  Pore  with  both  ice  and  liquid, 

pores. 


variations  in  pressure  and  geometry.  As  an  alternative 
to  this  approach,  we  can  assume  constant  vapor  pres¬ 
sure  and  look  at  the  dependence  of  equilibrium  tem¬ 
perature  on  variations  in  geometry.  When  all  three 
phases  of  water  are  present,  there  is  one  equation  (eq 
21 )  for  each  phase.  When  three  phase  boundaries 
occur,  there  is  one  equation  (eq  29)  for  each  boundary. 
These  equations  can  be  combined  in  various  ways  to 
derive  equations  which  are  applicable  to  particular 
situations.  For  example,  when  each  pore  in  a  soil  is 
either  i ce- or  water-filled  (Fig.  4), 


Jo  /  VJ° it  _  ^°g 

L  \  /"gs  '"eg 


(58) 


Thus  the  temperature  at  which  the  system  comes  to 
equilibrium  is  determined  by  the  radii  of  both  the  ice/ 
vapor  and  liquid/vapor  interfaces  in  the  pores. 

In  a  pore  containing  both  solid  and  liquid  (Fig.  5), 
the  equilibrium  temperature  and  geometry  are  related 
bv 


(59) 


where  Lt  is  the  latent  heat  of  evaporation. 

These  are  onlv  two  of  many  equations  which  can 
be  derived  tor  a  variety  of  situations.  Note  that  to 
define  equilibrium  unambiguously,  it  is  necessary  to 
specify  the  nature  of  the  interfaces  actually  existing 
in  the  soil. 

Grain  growth  in  snow 

In  snow  at  a  subfreezing  temperature,  the  vapor 
pressure  over  an  ice  surface  depends  on  the  radius  of 
curvature  of  the  surface.  (In  the  absence  of  an  im¬ 
posed  temperature  gradient,  snow  is  nearly  isothermal, 
so  Kelvin’s  equation  applies.)  By  analogy  with  eq  42, 
the  vapor  pressure  over  the  curved  surface  of  a  snow¬ 
flake  is 


ln(pg/p0) 


RT 


(60) 


so  larger  vapor  pressures  prevail  over  surfaces  of 
smaller  particles.  As  the  water  vapor  is  removed  by 
condensation  on  the  larger  particles,  the  smaller  parti¬ 
cles  evaporate;  grain  growth  proceeds  by  this  mechanism. 

When  the  pore  space  of  snow  is  filled  with  liquid 
rather  than  vapor,  the  fluid  pressure  is  essentially  con¬ 
stant  but  large  temperature  differences  arise  among 
the  particles.  For  a  constant  liquid  pressure,  the  tem¬ 
perature  of  each  particle  is  given  by 


27q 

U>,  rp 


(61) 


where  rp  is  the  radius  of  the  ice  particle.  Since  the 
smaller  particles  are  at  lower  temperatures  than  the 
larger  particles,  heat  flows  .owards  the  smaller  parti¬ 
cles,  causing  them  to  melt,  while  the  larger  particles 
grow  by  accretion.  This  mechanism  causes  rapid  grain 
growth  and  grain  rounding  in  liquid-saturated  snow. 


Liquid-filled  veins 

Metallurgists  discovered  some  time  ago  that  liquid- 
filled  veins  can  occur  at  three-grain  boundaries  in 
crystalline  materials  (Fig.  6).  These  veins  also  occur  in 
ice,  whose  dihedral  angle  is  about  20°.  Veins  provide 
a  fluid  path  and  render  ice  permeable,  even  at  subzero 


Figure  6.  Liquid -filled  vein  at  a  three- 
groin  boundary. 
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Figure  7.  Equilibrium  temperature  of  a  liquid 'filled 
vein  as  a  function  of  vein  size. 


Figure  8.  Ice  particle  enclosed  in  liquid  (a) 
and  a  separate  ice  particle  and  water  droplet  (b). 


temperatures.  The  temperature  of  a  two-dimensional 
vein  is  given  by 

Tm  =-0.137  (62) 

p^L  r v 

where  the  coefficient  is  determined  by  the  geometry 
and  rv  is  the  radius  of  the  largest  circle  which  will  fit 
in  the  vein.  The  temperature-radius  relation  shown  in 
Figure  7  is  an  important  part  of  the  “candling”  of 
lake  ice.  As  the  internal  ice  temperature  rises  by  radi¬ 
ation  absorption  and  heat  conduction,  the  veins  en¬ 
large.  Thus,  grain  boundary  strength  decreases  and 
permeability  increases  with  increasing  temperature. 

Ice  clouds 

As  stated  earlier,  equilibrium  can  only  be  defined 
unambiguously  when  the  nature  of  the  interfaces  is 
known.  For  example,  ice  particles  in  clouds  can  grow 
either  from  the  vapor  or  the  liquid  (Fig.  8).  When  an 
ice  particle  is  surrounded  by  liquid,  no  solid/vapor  in¬ 
terface  exists,  so  one  less  equation  is  available  to 
specify  equilibrium.  In  this  case  the  three  equations 
from  eq  21  and  two  equations  from  eq  29  are  com¬ 
bined  to  give  the  temperature  at  which  equilibrium 
occurs,  yielding 

Tm-  jte)qii-^Jr0.  (63) 

L  I  rtt  rii  \ 


When  the  liquid  droplets  and  ice  particles  are  separate, 
the  three  equations  from  eq  21  and  three  equations 
from  eq  29  are  combined  to  give  the  equilibrium 
temperature 


=  2  f^ggg 
L  rCg 


^sg 

rsg 


To¬ 


rn 


These  two  cases  differ  significantly,  since  the  effect  of 
the  radius  of  the  liquid  droplet  is  very  small  (and  nega¬ 
tive)  in  the  first  case  but  large  (and  positive)  in  the 
second  case. 


CONCLUSION 

The  principles  given  here  can  be  applied  to  a  wide 
range  of  cold  regions  problems.  These  basic  ideas  con¬ 
cerning  the  laws  of  thermodynamics  and  Laplace’s  equa¬ 
tion  explain  a  large  number  of  capillary  phenomena.  No 
doubt  many  future  observations  will  be  enriched  by 
these  well  tested  concepts. 
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